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We discuss the computation of transport coefficients in large Nf QCD and the 
1 0(N) model for massive particles. The calculation is organized using the 1/7V 

nJ ' expansion of the 2PI effective action to next-to-leading order. For the gauge theory, 

^_>^ ' we verify gauge fixing independence and consistency with the Ward identity. In 

the gauge theory, we find a nontrivial dependence on the fermion mass. 



1. Motivation 

In understanding the evolution of quantum fields out of equilibrium, ther- 
malization plays an important role. For long enough times, the dynamics of 
a typical system can be characterized by transport coefficients, as they con- 
$H ' trol the relaxation of a system towards equilibrium. For example, in heavy 

^ ion collisions, the presence of an appreciable shear viscosity would modify 

the ideal hydrodynamical description^ of the aftermath of the collision. ^'^ 
The nonequilibrium evolution of quantum fields and the subsequent 
thermalization have been studied using various approximations: mean-field 
dynamics (such as Hartree or leading-order large iV approximations) , inho- 
mogeneous mean-field theory, the classical approximation, kinetic theory, 
dynamics based on truncations of the 2PI or 2PPI effective action, etc. It is 
clear that any scheme aimed at describing the long-time behavior should, 
among other things, describe correctly the approach to equilibrium of a 
typical system. One requirement for an approximation scheme is there- 
fore to yield sensible results for transport coefficients. While it is clear 
that some approaches (mean-field theory, classical approximation) do not 



* Based on presentations at Strong and Electroweak Matter (SEWM04), Helsinki, Fin- 
land, June 16-19 2004, and the Workshop on QCD in Extreme Environments, HEP 
Division, Argonne National Laboratory, USA, June 29-July 3 2004. 
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give the correct long-time quantum behavior, for others this might be more 
involved. 

One approach to the dynamics of quantum fields far from equilibrium 
that has been successful in the past few years*" employs the 2PI effective 
action. It is therefore interesting to see whether the 2PI framework gives 
a reliable description of transport coefficients, in the limit where a semi- 
analytical computation can be carried out, i.e. in a weak coupling or a 
large N expansion. Below we review how transport coefficients can be 
computed using the 2PI effective action as an organizational tool.^^ We 
discuss explicitly two theories in the large N limit: the shear viscosity in 
the 0{N) model^^ and the electrical conductivity and shear viscosity in 
large Nf QCD with massive fermions.^® 

2. 2PI effective action 

The 2PI effective action offers an approach to deal with the infinite hierar- 
chy of correlation functions in field theory, based on a variational principle 
for one- and two-point functions. For a bosonic field, with vanishing mean 
field {(j)) = 0, the effective action is written as^^ 

T[G] = |TrlnG-i + '-Tr G^'iG - Go) + T^[G], (1) 

where Gq^ is the free inverse propagator. T2[G] contains all two-particle 
irreducible diagrams without external legs. For fermionic fields the I's 
arc replaced by — I's. Although the attention is usually focused on the 
(one- and) two-point functions, the 2PI effective action is formulated as 
a generating functional with external sources, similar to the IPI effective 
action, and therefore higher order n-point functions are accessible as well. 
In particular, there is a four-point vertex function obeying a Bethe-Salpeter 
equation, which in momentum space roads 

r(*)(P, K) = A(P, K) + ]- [ A(P, R)G^{R)r'-''\R, K). (2) 
^ Jr 

This equation sums ladder diagrams with a kernel A obtained from the sec- 
ond derivative of r2 [G] . It plays therefore a crucial role in the computation 
of transport coefficients^^ in the case that those arc dominated by ladder 
diagrams, such as in scalar^"'-^* and largo Nj gauge theories^^ to leading 
order and in weakly coupled gauge theories to leading logarithmic order. 

We now apply this to two theories in the large N limit: the 0{N) model 
and large Nf QCD. We consider the shear viscosity in the 0{N) model with 
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massive particles both at weak coupling as well as in the large N limit. 
For a weakly- coupled single-component scalar field the shear viscosity was 
computed some time ago by Jeon.^° In large Nf QCD, we present the first 
diagrammatic calculation of the electrical conductivity and shear viscosity 
for massive fcrmions. For massless fcrmions, these transport coefficients 
were computed before by Moore^^ using kinetic theory. 




Figure 1. 2PI contribution to the cfTcctivc action at LO and NLO in tiie large N limit 
in the 0{N) niodi-l (k-Il) and at NLO in large Nf QCD (right). The dashed line in the 
0(N) model sums the chain of bubbles, see Fig. 2. 

The starting point is the 2PI contribution to the effective action in 
the 1/N expansion, shown in Fig. 1. In the 0{N) model there is a LO 
and a NLO contribution.^'^ In large Nj QCD there is only a contribution 
at NLO. The chain of bubbles in the 0{N) model is indicated with the 
dashed line (see Fig. 2). In this formulation the scalar and gauge theory 
arc conveniently similar. Wo emphasize that all results below arc generated 
from the graphs in Fig. 1, indicating the power of the 2PI formalism. 

>-< = X*X)-< 

Figure 2. Auxiliary correlator in the 0{N) model, summing the chain of bubbles. 

Extremizing the effective action yields the Schwinger-Dyson equations 
for the two-point functions. For the scalar case we find 

G-i=Go-i-S, D-'=Do'-Il, (3) 

where G is the scalar and D the auxiliary correlator. For the gauge theory 

5-i = 5o-i-I], D-' = D^'-n, (4) 

with S the fermion and D the gauge field propagator. The self energies, 
depending on full propagators, are shown in Fig. 3. In the computation of 
transport coefficients, it is crucial to use dressed propagators for two rea- 
sons: to screen the so-called pinching poles, ^"^'^^ reflecting the dependence 
on the finite lifetime of quasiparticles due to collisions in the plasma, and 
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Figure 3. Self energies at LO and NLO in the large N limit in the 0{N) model (left) 
and at NLO in large Nf QCD (right). 

to screen the divergences due to the exchange of offshell gauge bosons^'^ 
(this only in the gauge theory). 

Differentiating the self energies yields the rungs that appear in the 
Bethe-Salpeter equation. We find that some rungs contribute to the trans- 
port coefficients at leading order, whereas other rungs are subleading and 
can be neglected. When considering the Bethe-Salpeter equation for other 
purposes (e.g. for renornialization^''''^^'^^) wc expect that all rungs that fol- 
low from the effective action might have to be considered. 
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Figure 4. Kernels at LO and NLO in the large N limit in the 0{N) model (left) and 
at NLO in large Nf QCD (right). 



In order to distinguish leading from subleading rungs, we need to dis- 
cuss power counting. In the large TV expansion, this is fairly straightfor- 
ward. Positive powers of N arise from two sources: closed loops of scalar or 
fermionic lines and scalar or fermionic propagators suffering from pinching 
poles (pinching poles are screened by the scalar or fermionic thermal width 
Fp and result in contributions ~ l/Fp ~ Negative powers of N arise 
from the coupling constants, which are taken to scale as 1/A'' as is usually 
done in 1/A^ expansions. As a result we find the rungs given in Fig. 4. In 
the 0{N) model, the first point-like rung does not contribute for kinemati- 
cal reasons. In large Nf QCD, the contribution from generic onshell gauge 
fields is subleading: the gauge field propagator appears therefore only in 
rungs and there is no need to consider the thermal width of onshell gauge 
bosons. The thermal width of an onshell fermion determined by the self 
energy in Fig. 3 is ill-defined. ^^'^^ However, the problematic part cancels 
against part of the contribution from the line diagram in the kernel. ^'^'^^ 
Finally, we note that the two box diagrams in large Nf QCD differ only 
in the orientation of the fermion lines; this ensures that Furry's theorem is 
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satisfied. 



3. Summation of ladders 

To sum the ladder series, we use the technique recently presented by Valle 
Basagoiti,^^ employing the Matsubara formalism and an efi'ective three- 
point vertex T>. In terms of this vertex, the integral equations are shown 
in Figs. 5, 6. Note that the propagators in these diagrams are still dressed. 




Figure 5. Integral equation for tlie effective vertex in the 0{N) model. The first closed 
loop on the RHS does not contribute. 



We consider these equations in the kinematical configuration special for 
transport coefficients: the momentum entering on the left is Q = (g" — > 
0, 0), while the momentum entering and leaving on the right is P = (wp, p), 
where ojp = -^/p^ + rT? with m the scalar or fermion mass. Transport 
coefficients are then extracted from the correlator obtained by closing the 
lines with an insertion of the appropriate current, indicated with the small 

black dot.21.17,18,19 




■ 00000 




. 00000 . 






Figure 6. Integral equation for the effective vertex in large Nf QCD. In the case of the 
electrical conductivity, the last two contributions cancel. 



The basic quantity we use to solve the integral equations is the ratio of 
the effective vertex V{p) and the thermal width Pp. For the shear viscosity 
we define 

Wp Fp [Cf/Tf large iV/ QCD 



with Cf = {N^ ~ l)/2iVc, Tf = 1/2. Introducing the color factors in 
this way removes all color factors from the integral equation. The shear 
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dpTp 2Nf [°° p 



'^'-^Y^i "o^/'Mxip) (6) 

with dp = In terms of % the integral equations read compactly as 

f-OC 

Hp)x{p) = S{p) + / drn{p, r)xir), (7) 



where J^{p) oc Fp, S{p) is determined by the bare vertex, and 7i — Tiuno + 
Hbox is determined by the rungs. Since 7i{p,r) = T-l{r,p), the integral 
equations follow as the extremum condition from the frmctional 

/■oo r ^ I f°° 

Q[x]= dp S{p)x{p) - -^:F{p)x^{p) + - drH{p,r)x{r)x{p) , 

(8) 

which allows for a variational treatment. The value of Q at the extremum 
is immediately proportional to the transport coefficient. 



4. Gauge fixing and Ward identity 

In applications of 2PI effective action techniques to gauge theories, gauge 
fixing and Ward identities have to be considered. ^^^■^"^'^^ It is therefore in- 
teresting to see where gauge fixing parameters appear and why they drop 
out in the calculation. We use the generalized Coulomb gauge such that 
the gauge field propagator, 

pH pv 

= ATPr + ^i5''V"+^^^, (9) 

consists of a transverse, a longitudinal, and a gauge fixing piece. The gauge 
fixing parameter ^ appears in three places: 

(1) thermal width. The imaginary part of the fermionic self energy 
is proportional to the discontinuity (spectral density) of the gauge 
boson propagator, i.e. to pr, Pl- Since the gauge fixing part has no 
discontinuity, ^ drops out. 

(2) line diagram. In the kinematical limit we consider also this con- 
tribution is proportional to the discontinuity of the gauge boson 
propagator and ^ drops out. 

(3) box diagram. In the kinematical limit we consider all fermion lines 
are onshell and ^ drops out exactly. 
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One can also verify that the effective vertex in Fig. 6 and the fermion 
self energy are related via the Ward identity, again in the kinematical limit 
we consider. The analysis in this case is in fact easier than for the weakly 
coupled case,^^ since the contribution to the thermal width arising from 
soft fermions in the latter is now subleading in the large Nf expansion. 

When extending this approach to dynamics far from equilibrium it is 
necessary to include contributions that could be dropped in the leading 
order analysis of transport coefficients presented here. The gauge fixing 
dependence that will be present in that case should be suppressed for suf- 
ficiently large Nf. 



5. Variational solution 

The integral equations are in general too complicated to be solved analyti- 
cally. In the scalar case however, we found an approximate but surprisingly 
accurate solution for zero mass and vanishing coupling in the limit of ul- 
trahard momentum p T. We find 

;imxW =6912.^^1^. (10) 

Using this result yields for the shear viscosity in the weakly coupled 0{N) 
model, 

27648C(5) iV3 t3 ^^^^ ^ 

In order to obtain the full N dependence and not just the leading order 
behavior in the large N limit, we used the three-loop expansion of 

the 2PI effective action in the 0{N) model. -'^^ The result for the numerical 
prefactor is very close to the full results obtained numerically by Jeon^° 
(3040) and Arnold, Moore and Yaffe^^ (3033.5) for TV = 1. 

For arbitrary values of the mass and coupling constant (limited by the 
presence of the Landau pole) , we solve the integral equations variationally. 
The function x is expanded in a set of trial functions and the remaining 
integrals (a one-dimensional integral for S, three-dimensional integrals for 
J-" and Tiiino, and a four-dimensional integral for Hhox) are performed nu- 
merically. We found that a set of four trial functions suffices. 

The shear viscosity in the 0{N) model, normalized with the analyti- 
cal result (11) in the large A'' limit, is shown in Fig. 7 as a function of 
renormalizcd mass at zero temperature for various values of the running 
coupling constant (the result is renormalization group independent). We 
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Figure 7. Shear viscosity in the 0{N) model vs. the renormahzed mass at zero temper- 
ature for various values of the coupling constant A(/i = T). 

find that the shear viscosity increases monotonically with increasing mass. 
This behavior was also found by Jeon^'' for A'^ = 1. Interestingly, we find 
surprisingly little dependence (apart from that contained in 7700) on the 
coupling constant, for all values allowed by the Landau pole. 

For large Nf QCD the shear viscosity is shown as a function of ferniion 
mass in Fig. 8, again for several values of the effective coupling constant 
5eff (Seff = '^pg'^Nf). The viscosity is normahzed with 770 = dp / (TfCf) x 
T^/g^. In this case we observe a nontrivial dependence on the mass. After 
a slight increase, we find that the viscosity decreases with increasing mass. 
This behavior is caused by longitudinal gauge bosons below the lightcone. 
In the limit of very large mass it can be shown'^'* that the viscosity goes to 
zero.*^ We also find a much stronger dependence on the coupling constant, 
compared to the scalar theory. 

6. Conclusion 

The diagrammatic calculation of transport coefficients in the 0{N) model 
and large Nf QCD is well organized when the 2PI effective action is used: 

'^This is not in contradiction with the conjectured lower bound on the viscosity/entropy 
ratio, '^^ since the entropy goes to zero exponentially. 
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Figure 8. Shear viscosity in large Nf QCD vs. the fermion mass for various values of 
the effective coupling constant gcfiifJ- = Trcxp"'''^ T). 

all necessary summations are generated by the few diagrams in Fig. 1. In 
the scalar theory the shear viscosity increases monotonically with mass and 
has a dominant 1/A^ dependence on the coupling constant. In the gauge 
theory, on the other hand, the dependence on the coupling constant and the 
fermion mass is highly nontrivial, due to longitudinal gauge bosons below 
the lightcone. 

We emphasize that understanding transport coefficients diagrammati- 
cally provides necessary insight in the dynamics of quantum fields out of 
equilibrium. Our results provide further support for the applicability of 
truncations of the 2PI effective action to noncquilibrium QFT. 
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